Assessing Alternatives for Directional Detection of a WIMP Halo 
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The future of direct terrestrial WIMP detection lies on two fronts: new, much larger low back- 
ground detectors sensitive to energy deposition, and detectors with directional sensitivity. The 
former can probe a large range of WIMP parameter space using well tested technology while the 
latter may be necessary if one is to disentangle particle physics parameters from astrophysical halo 
parameters. Because directional detectors will be quite difficult to construct it is worthwhile explor- 
ing in advance generally which experimental features will yield the greatest benefits at the lowest 
costs. We examine the sensitivity of directional detectors with varying angular tracking resolution 
with and without the ability to distinguish forward versus backward recoils, and compare these to 
the sensitivity of a detector where the track is projected onto a two-dimensional plane. The latter 
detector regardless of where it is placed on the Earth, can be oriented to produce a significantly 
better discrimination signal than a 3D detector without this capability, and with sensitivity within 
a factor of 2 of a full 3D tracking detector. Required event rates to distinguish signals from back- 
grounds for a simple isothermal halo range from the low teens in the best case to many thousands 
in the worst. 



> 
OS 

00 

o 
in 
o 

Or 
i 1 

o ■ 
a: 

CZ3 . 



INTRODUCTION 

Direct detection experiments for weakly interacting 
massive particles (WIMPs) continue to set ever more 
stringent limits on the nucleon-WIMP cross section 0, 
A precise understanding of the backgrounds is 
required to identify an excess of nuclear recoils. There is 
no unique signature that can separate neutron induced 
recoils from WIMP induced recoils in these detectors. 
Even annual modulation, which is at best a few percent 
effect, might be accounted for by seasonal background 
variations. (This is for a pure WIMP signal. When a uni- 
form background is included the effect is even smaller.) 
It has been recognized that a stronger signal comes from 
measuring the direction of the recoiling nucleus [6j al- 
lowing for a WIMP signal to be identified from a few 
events even in the presence of backgrounds 0, H| • Detec- 
tors that might measure the recoil direction have been 
designed and built. For example, time projection cham- 
bers have been used by DRIFT and NEWAGE 
and a scintillator with direction dependent response has 
also been studied 

In their present state the directional detectors are 
rudamentary at best. The current version of DRIFT, 
for example, does not measure the full three dimensional 
track of the recoiling nucleus. Instead it measures the 
recoil track projected onto a plane and has a number of 
other limitations. A detailed statistical study of DRIFT'S 
current capabilities has been performed |12| 

A new generation of directionally sensitive detectors 
are now being envisaged. Building a full three dimen- 
sional detector is a challenging, costly proposition. Is it 
necessary? Given the technical challenge and cost con- 



straints of these detectors, which changes will lead to the 
most sensitive detector? Here we provide a general anal- 
ysis of various design goals to determine the number of 
events required for detector designs ranging from two di- 
mensional to fully three dimensional detectors. We do 
not focus on any particular detector technology. Instead 
we apply a consistent set of parameters to a variety of 
detector configurations. This allows us to determine the 
optimal design goals. In the process we also are able to 
dispel various myths about optimal locations on Earth 
for detector placement, and the relative efficacy of two- 
dimensional versus three-dimensional event reconstruc- 
tion. 



THEORETICAL MODEL 



Detector Characteristics 

To quantify the capabilities of directionally sensitive 
detectors we consider a consistent, generic set of param- 
eters. For the detector target we use a xenon (mjsr = 
131 GeV) nucleus. We assume a threshold of Q t h = 
10 keV and two different WIMP masses, m x = 100 GeV 
(m x fts ton) and m x — 1000 GeV (m x 3> tun). For our 
Galaxy we focus solely on an isothermal model for the 
WIMP halo distribution, 



f(v) 



(1) 



We study three different values for vo spanning the 
range of current expectations, 170 km/s, 220 km/s, and 
270 km/s. The escape velocity of WIMPs from the 
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Galaxy is taken to be 650 km/s. The Earth's rotation 
axis is oriented at an angle e = 42° with respect to the 
Sun's motion. This value is relevant for the two dimen- 
sional detector We stress that these choices have been 
made to provide a consistent set of parameters to allow 
the intercomparison of detector designs. Thus, what will 
be important about our results will not be absolute con- 
straints, but relative ones, although we expect the overall 
order of magnitude for the required event rates will not 
differ compared to our estimates. 

Differential Event Rates 

The technique for calculating the WIMP scattering 
rate is well known 0- The differential rate as a func- 
tion of nuclear recoil direction (7, ip) is given by 

dR = ^oPo_ r d 3zj vF 2 (Q{Jv))f( z + z m y (2) 

Here mjy is the mass of the target nucleus, m x is the 
WIMP mass, v® is the velocity of the Earth through the 
WIMP halo, Q is the recoil energy of the nucleus, ero 
is the cross section for WIMP scattering off the target 
nucleus, and p is the local WIMP halo density. We 
consider only spin independent interactions and use the 
standard Helm form factor [HI for F 2 (Q). The geometry 
of the WIMP scattering gives 

Jv = v x sin 7 cos 4> + v y sin 7 sin <j> + v z cos 7. (3) 

The integration region, 1Z, is defined by the detector 
threshold, Qth, at the lower limit, 

y 2m x m N 

and the galactic escape velocity, v csc , at the upper limit, 

K + %,,) 2 + (v y + v^^) 2 + {v z + v @ . z ) 2 < v 2 csc . (5) 

See for a more detailed discussion. A full three dimen- 
sional detector probes the reaction rate outlined here J5J). 

Measuring the direction the nucleus is moving along a 
track is not always possible. For a three dimensional de- 
tector without forward-backward discrimination, a recoil 
in a direction (cos 7, </>) cannot be distinguished from a 
recoil in a direction (— cos 7, + 7r). The event rate for 
these two directions thus combine giving a total rate for 
the direction (cos 7, 4>) of 4^ (cos 7, (f>) + 4§ (— cos 7, 0+7r) 
where the differential rates are again given above (J2J). 

A two-dimensional detector can only resolve the recoil 
direction projected onto a plane. We assume for com- 
parison purposes however that it can be designed with 
forward-backward discrimination. Suppose the normal 



to the detector plane is oriented at angles (£, rj) with re- 
spect to the direction of the Sun's motion. The rate is a 
function of a single angle 4>' measured in this frame, 

dR = ^0A L _ f X d{cQsi) f d 3v> { J v y F 2 { Q< )f{ y + v e ). 

d<p' Trm N m x J_x J n 

(6) 

Here, primed coordinates are measured in the frame of 
the detector (where the normal vector points along the 
z-axis) and unprimed coordinates are measured in the 
frame of the Sun's motion. These two frames are related 
by rotations through the angles £ and r\. 

For a detector fixed to the Earth's surface, the detec- 
tor orientation with respect to the Sun's motion changes 
throughout the day due to the Earth's rotation. Let 9 be 
the angle between the detector's normal and the Earth's 
rotation axis and e be the angle between the rotation axis 
and the Sun's velocity. Then 

cos Q — cos e cos 9 + sin e sin 9 cos(27rf) (7) 

for t measured in days. 

It is also important to recognize that even full three- 
dimensional detectors will not have perfect angular reso- 
lution. To model realistic angular resolution either due to 
dispersion of the recoiling nucleus along its ideal track or 
due to inherent precision of the inherent detection mech- 
anism itself we convolve the ideal scattering rate @ with 
a smoothing kernel K (0, 17'), 

d ^m-Jm^(n'w. ( 8 ) 

We use a Gaussian smoothing about the direction of the 
ideal recoil, 

K(fl, n') = e-« 2 / 2CT V(2^) 3/ Verf(\/2 ( 7), (9) 

where 

cos£ = sin 9 sin 9' cos(</> — <fi') + cos 9 cos 9' . (10) 
We study this as a function of the width of the Gaussian. 

RESULTS 

We test the capabilities of each type of WIMP detec- 
tor by assessing their ability to distinguish the WIMP 
distribution from a flat background, using an isothermal 
halo as a fiducial test model. The probability that a 
WIMP will recoil in a particular direction, tti, is given 
by P(S!;) = ^f(^;)j where we have normalized the rate 
such that R= 1. Thus, we are probing the shape of the 
recoil spectrum. The likelihood function for N e detected 

events is defined by C = P(S1,:). We generate at least 

i=i 
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TABLE I: The number of events required to identify a WIMP 
signal above a flat background for different types of detectors 



and a WIMP mass of m v 



100 GeV. 



Detector 




v (km/s) 




Type 


170 


220 


270 


3D (full) 


6 


11 


18 


3D without FB 


176 


1795 


> 35,000 


2D — best / worst 


19/45 


34/75 


61/123 


2D rotating 


13 


24 


43 




TABLE II: The number of events required to identify a WIMP 
signal above a flat background for different types of detectors 
and a WIMP mass of m x = 1000 GeV. 



Detector 




v (km/s) 




Type 


170 


220 


270 


3D (full) 


14 


27 


51 


3D without FB 


152 


217 


371 


2D fixed — best/worst 


51/129 


97/217 


175/368 


2D rotating 


31 


61 


125 



100,000 sample distributions for each N e and apply the 
log-likelihood test to find the minimum number of WIMP 
events such that we have a 95% detection 95% of the time 
(see @ for more details). 

The results for m x — 100 GeV are given in Tabled for 
the range of detectors we have considered, where 'full' re- 
flects a full three-dimensional detector with perfect angu- 
lar resolution. We shall discuss the degradation implied 
by limited resolution shortly. The same set of results for 
m x — 1000 GeV are given in Tabled Although we have 
restricted our quantitative study to isothermal models 
the qualitative features of the comparison remain valid 
for other models, including models dominated by single 
streams of WIMPs . 

Our results underscore the need for forward-backward 
detection. Indeed, this is the single most important fea- 
ture that allows directional detectors to gain sensitivity 
to the WIMP signal compared to backgrounds. Since 
spin independent WIMP scattering is azimuthally sym- 
metric about the direction of the incoming WIMP the 
dominant WIMP signal comes from the a comparison of 
forward-backward scattering events. This is seen in the 
results in Tables |U and [n] A three dimensional detec- 
tor, even with perfect angular resolution, but without 
forward-backward discrimination requires at least an or- 
der of magnitude more events than a three dimensional 
detector with such discrimination and even many more 
than a poorly aligned two dimensional detector to distin- 
guish a WIMP signal from terrestrial backgrounds. This 
is because without forward-backward discrimination the 
detector relies upon the difference between head-on and 
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FIG. 1: The number of events required as a function of 
the full- width half maximum (FWHM) of the smoothing ker- 
nel @ for the isothermal models and m x — 100 GeV. 
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FIG. 2: The number of events required as a function of the 
angle between the detector normal and the direction of the 
WIMP wind (Sun's motion), 9 for the isothermal models and 
m x = 100 GeV. 

glancing (wide angle) collisions as well as high angular 
resolution to distinguish a WIMP signal from the back- 
ground. 

We next explore how the sensitivity of a three- 
dimensional detector depends upon its angular resolu- 
tion. In figure ^ we display the number of events re- 
quired as a function of the angle of the full-width half 
maximum (FWHM) of the detection cone for the events. 
Note that as the angular resolution degrades, the num- 
ber of events required for a three dimensional detector 
quickly approaches that of a two-dimensional detector, 
as expected. In order to be significantly more efficient, 
the angular resolution of such a detector must be better 
than about 60 degrees (fwhm). 

We finally focus on two-dimensional detectors, in part 
because these are likely to be the most practical in the 
near future, and because less attention in the past has 
been paid to them than hypothetical three dimensional 
detectors. 
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It is clear that the efficacy of a planar detector will 
depend upon the orientation of its plane with respect 
to the direction of the WIMP wind. Specifically a two- 
dimensional detector fixed to the Earth will be oriented 
so that its normal vector makes an angle 9 with the 
Earth's rotation axis. The choice of 9 determines how 
much time the detector will spend at various angles C, 
relative to the Earth's direction of motion. An orienta- 
tion of C = 0° is clearly the worst since the detector plane 
is then perpendicular to the WIMP wind. The number 
of events required is a function of the angle 9 chosen for 
the detector as shown in figure |2] 

Tables [I] and [H] give the minimum and maximum num- 
ber of events required for optimal versus worst-case ori- 
entation of the detector. Note that the shape of the func- 
tion in figure |21 depends on the orientation of the Earth's 
axis, e = 42° relative to the motion of the Sun through 
an isotropic halo. For halos in which a WIMP stream 
arose which was not due to the motion of the Sun, the 
shape of the curve would change. Low values of e in this 
case would produce a more steeply descending function, 
and high values would produce an ascending function. 

Note that the results thus far assume that the detec- 
tor axis is fixed to the plane of the Earth. In this case 
the time averaged number of events required is on aver- 
age 3 times less for a two-dimensional detector oriented 
in the best possible axis (ie. at or near 9 = 90° for 
isothermal halo distribution). Such a detector however 
requires three to four times the number of events of a full 
three-dimensional detector. Whether or not achieving 
this additional factor of 3 in sensitivity for a full three- 
dimensional detector is worth the technological challenge 
is not clear. However, as seen in the tables, this fac- 
tor of three can be reduced by using a detector which is 
not fixed relative to the plane of the Earth, but which 
can rotate over the course of each day with respect to 
the earth to maintain an optimal orientation with re- 
spect to an expected WIMP wind. For an isothermal 
halo, roughly only 2 times as many events are required 
for such a rotating detector compared to a detector with 
full three-dimensional tracking capability. The technical 
difficult of producing the former may be less demanding 
than that required to produce the latter. 

Our results can be summarized as follows: For di- 
rectional WIMP detectors, forward-backward discrimi- 
nation is far more valuable than three-dimensional res- 
olution of the track, at least for the isothermal model 
considered here. Furthermore, a two-dimensional detec- 
tor can be oriented into the predominant WIMP wind 
so that the number of events required to distinguish a 



WIMP halo from a terrestrial background is comparable 
to that required for even a full three dimensional detec- 
tor. While there may be other reasons (including back- 
ground rejection) for considering fully three dimensional 
directional detection methods, our results thus suggest 
that concentrating simply on forward-background sensi- 
tivity is the most important new direction that should 
be pursued in directional WIMP detection, and that pla- 
nar detectors can, in this case, provide nearly optimal 
detectional sensitivity (We note that as this work was 
being prepared for submission a related study appeared 
which complements this work, and is consistent with our 
results 
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